Low temperature dipolar echo in amorphous dielectrics: Significance of relaxation and 

decoherence free two level systems. 
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04 \ The nature of dielectric echoes in amorphous solids at low temperatures is investigated. It is shown 

that at long delay times the echo amplitude is determined by a small subset of two level systems 
r-{ ' (TLS) having negligible relaxation and decoherence because of their weak coupling to phonons. The 

echo decay can then be described approximately by power law time dependencies with different 
powers at times longer and shorter than the typical TLS relaxation time. The theory is applied to 

■ recent measurements of two and three pulse dipolar echo in borosilicate glass BK7 and provides a 
perfect data fit in the broad time and temperature ranges under the assumption that there exist 
two TLS relaxation mechanisms due to TLS-phonons and TLS-TLS interaction. This interpretation 

, - is consistent with the previous experimental and theoretical investigations. Further experiments 

^ ■ verifying the theory predictions are suggested. 

' ^ . Two level systems (TLS) in low temperature amorphous solids have recently attracted growing attention due to 
their performance limiting effects in superconducting qubits for quantum computing^""^ and kinetic inductance photon 
detectors for astronomy*. TLS are represented by atoms or groups of atoms tunneling between two close energy 
minima (see Fig. [H Refs.^'^). Understanding of the TLS effect on thermodynamics and kinetics of materials requires 
knowledge of their dynamical properties including relaxation and decoherence rates. These rates can be determined 
' using the spin-echo technique^ generalized for amorphous solids.^ However, the time-dependence of echo signals is 
51^ ] rather complex and cannot be expressed using the simple exponential decay directly determining TLS relaxation or 
O ■ decoherence rates. Instead, since all TLS interact differently with the external pulse the decay of echo signal as a 
, ^ , ' function of delay time is much more complicated 4"— To our knowledge there is no general theory capable to describe 
the complicated echo signal and use it to extract the information about TLS relaxation and decoherence. Below we 
T-H ■ propose such a theory averaging individual TLS responses over their parameters with the special attention to their 
^ \ coupling to phonons. In our scenario the echo signal is determined by the small subset of relaxation and decoherence 

■ free TLS formed due to fluctuations of environment in a certain extent similarly to Refi^i, where this idea has been 
] proposed and explored for quantum bits (qubits). A common sense based assumption about TLS coupling statistics 

■ results in an almost perfect fit of experimental data for Borosilicate glass BK7 with the fit parameters consistent 
with the previous experimental data obtainedJ^ In these experiments the echo amplitude has been observed for 

. unprecedentedly long times following its five order of magnitude decay (see Figs. [3] and H]) which creates necessary 

■ grounds for the present theoretical development. 

I , TLS can be characterized by a potential well asymmetry A and a tunneling amplitude Aq, distributed according to 
• • ■ the universal law, P(A, Ag) ~ Pq/Aq, which reflects the exponential sensitivity of a tunneling amplitude to a potential 
barrier. TLS is also coupled to continuum environment including phonons and other TLS. For each TLS, i, its coupling 
to environment is determined by its elastic tensor 7^ describing its interaction with the strain field (phonons) , e, given 
^ ' by —^f^eabSf {a,b = x,y,z). It is convenient to represent the TLS elastic tensor by its invariant | 7 |« yTr^, 



having average root mean squared value 70 because the limit I7I — >■ 0, where TLS relaxation and decoherence rates 
also becomes very small compared to their average values, is relevant for the long time echo amplitude as pointed out 
below. ^ 

For the further consideration one can break the TLS elastic tensor into two parts including the transverse, j\ = 
7 — iTr(7)/ and longitudinal, 7'; — /Tr(7), parts (/ is the 3x3 unit matrix). The longitudinal part is coupled only 
to longitudinal phonons having larger sound velocity, c;, than the transverse ones, Ct- Since the relative weight of the 
longitudinal contribution to TLS interaction and relaxation decreases with the parameter ct/c; (for instance for the 
phonon stimulated relaxation it scales as2^ (ct/c/)^ 1) one can approximately ignore the longitudinal part of the 
elastic tensor and restrict the consideration to its transverse part (see discussion after Eq. (jH])). 

TLS interaction with phonons and with each other results in their relaxation to equilibrium which is the subject 
of interest for echo experiments. According to the earlier theoreticaU^"— and experimentali^"— studies the TLS 
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relaxation rate in dielectric glasses can be represented as 
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where £' = -^Z + Aq is the TLS excitation energy, p is the material density, c is the characteristic speed of sound, 
and Uq is the average absolute value of TLS-TLS elastic interaction constant for 1/ interaction, while their dipole- 
dipole interaction is usually smaller and neglected herei^ii^i^ It is convenient to characterize the TLS-TLS interaction 
constant by its average absolute value, since it enters in that form to the spectral diffusion22, and relaxation— The 
rate ks describes the relaxation of a "typical" symmetric (|7i| = 70, A = 0) TLS. The first contribution to the TLS 

relaxation rate ('^) '^'^^^ ^ 2kBT ^ '^^'^ TLS-phonon interaction, while the second one {BT) is attributed 

to the TLS-TLS interaction. Both experiment and theory are qualitatively consistent with Eq. ([T]), but the observed 
numerical prefactor B is larger by around two orders of magnitude than the related theoretical predictions^ 

The TLS i interaction with other TLS's is determined by its interaction constant U^i which determines the spectral 
diffusion and decoherence for the given TLS i?^ In the case of interest where TLS coupling to phonons is small 
I 7i |— this interaction scales linearly with | 7^ | (therefore the TLS-TLS interaction stimulated relaxation rate 
scales as I 7i p in Eq. (HI), cf-^^) so we can set Uoi — C^o^- 




FIG. 1: Model of TLS. A is the energy difference between left and right well states when isolated, which are coupled with the 
tunneling amplitude Aq. 

The TLS contribution to the echo amplitude decreases exponentially with the characteristic times of experiment, 
ti2, ^13 (see Eq. ([TSl) below). However, a subset of TLS having very small coupling to phonons, I7I ^ 70, can yet 
strongly contribute to the echo at arbitrarily long times. To describe the response of such subset we introduce the 
distribution function F{x) for the relative TLS-phonon coupling strength x — |7|/7o- In our model the interaction 
tensor 7 is a symmetric 3x3 traceless matrix having 5 independent components. To have I7I approaching one needs 
all these five components to approach zero. This suggests F{x)dx cx dx^ for x -C 1. At large a; > 1 we expect that 
its statistics approximately obeys the Gaussian distribution in agreement with the law of large numbers for the sum 
of independent contributions. Combining two limiting regimes we propose the distribution function for the relative 
coupling constant x in the form 

, 25710 4 5.2 
F{x) = ^^x^e- — . 2 

This function is chosen to be normalized by unity and to have the average squared value of x to be also equal to unity. 

Echoes occur when a series of pulses at the same frequency, Wq, are applied to the glassy material. Here we consider 
dielectric echoes formed by microwave pulses characterized by electric field amplitude, F, interacting with TLS dipole 
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moments, p ~ ID. Two pulse echoes are composed by two pulses of duration t ("7r/2-pulse") and 2t ( "7r-pulse" ) , 
separated by the time ti2, while the three pulse echoes are composed by three 7r/2-pulses, separated by times ti2 and 
i23 (see Ref..^>22 and Fig. [2]). 
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FIG. 2: Two (left) and three (right) pulse echo sequences 



The echo signal is observed at the frequency of the pump pulses after the time ii2 from the last pulse. This 
signal is caused by the reversed precession of the resonant TLS polarization {Eq — ftwo) turned by the tt angle due 
to a TT pulse (two pulse echo) or two 7r/2- pulses (three pulse echo, see Refsj^"— ). Actual "turn angles" of resonant 
TLS wavefunction during the pulse, (f = TpFAo/{2Eh) differ for different TLS and they are not equal to 7r/2 or tt. 
We assume that all turn angles are smaller than 1 so that the perturbation theory with respect to the external field 
is applicable, though this is not significant for the echo signal at intermediate and long times. The evolution of TLS 
wavefunction (ci, C2) initially occupying its ground state ci = 1, C2 = (consideration of excited state simply adds the 

factor of the equilibrium population difference tanh ^ 2^t ) ) during the two pulse echo sequence can be approximately 

described as§^ (cf. Eq. ^) 
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where the phases, $, are defined as h^i{ti2) = /J^^ dTE{T), h^2{ti2) = Jq^^ dTE{T + ti2), E{t) is the TLS energy, 
fiuctuating with the time due to spontaneous transitions of neighboring TLS. The echo signal can be observed at the 
time t'^2 — ^12 due to the approximate compensation of phases $1 « $2- Its intensity is given by the average TLS 
dipole moment off-diagonal element, oscillating with the resonant frequency. 



< P2p >« P^{CIC2 + C^Ci) OC - Sin(c.i)(l - 52)2g-fc_(l-,^).^t,. ^ > . 



(4) 



We assumed that the absolute value of TLS dipole moment is approximately constant, p} As it follows from the 
preliminary analysis this assumption results in the best fit for the echo amplitude dependence on the pulse time. The 
discussion of this problem in detail is beyond the scope of this paper. 
Similarly one can obtain the three pulse echo amplitude in the formi^^ 



< P3p >' 



sin(tt't)(l 



s{l-q^)x^(tl2+t23) ^ gi(*l(tl2)-*3(tl2,t23)) > 



2n$3(ii2,i23) 



dTE{T+t23)- 
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Echo amplitudes decay with the delay times, ti2, ^23, due to TLS relaxation, Eq. ([T]) and phase decoherencC) 
occurring due to its phase fluctuations associated with its interaction with neighboring TLS making transitions between 
their ground and excited states. These fluctuations lead to the decay of neighbor averaged phase exponents in Eqs. 
O, © as 



< e 



< g.(*i(ti2)-*2(ti2)) >^ exp [-xqG2p{ti2)] 

^(*l(tl2)-*3(t.2,t23)) [-XqGsp{h2,t23)] 



(6) 



The factor q appears in the exponents because the TLS "diagonal" interaction with neighbors scales as [/ cx while 
the factor x is because the elastic interaction coupling constant is proportional to the absolute value of TLS-phonon 
interaction tensor (see Ref.^^'^^). The linear dependence of exponents in Eq. ^ of these factors is a consequence of 
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the specific of tlie TLS interactionjiSiS^i^ As it is sliown in Sec. |T]the functions G'2p.3p can be approximated by 

the following analytical expressions similar to the widths of the spectral diffusion^ 



G2p{tl2) 



G3p(il2, ^23) — 



afc^rtia In (1 + r^apfcy^iz) _ 

fl ??2p 

_ akBTti2 In (1 + rj^pkritu + ^23)) . 



V3p 

-^BT.a^—PoUo. 



(7) 



Here kx is the characteristic relaxation rate of symmetric TLS (Aq — E ^ kgT, I7I = 70), the free parameter a ex- 
presses the interaction strength of TLS and the free parameter ^ accounts for the difference between average relaxation 
rates caused by TLS-TLS and TLS-phonon interactions (see Appendix). At short times ti2, ^23 <C l/kr we obtain 
G « '^PoUokTti2{ti2+t23) ^ j^jj accord with Refs j^^i^^ . The logarithmic behavior at longer times is due to integrated 
contribution of neighboring TLS possessing small tunneling amplitudes, Aq ^ kgT. The parameters r]2p, rj^p are 
used as the free parameters of the theory, describing the transition between two asymptotic time dependence. These 
parameters are determined below by the direct comparison of theory and experiment. The relaxation parameters, A 
and B, Eqs. (P), ([3]), are determined similarly. 

The integrated echo signal is calculated by averaging of Eq. ([7]) over the distribution of TLS tunneling amplitudes 
and the statistics of their coupling constants, Eq. ([2]). The final result for the echo signal can be expressed as 



hp,3p oc / X e ^ dx dq{l - q )e 



»fcgrti2 '"(l + 12p.3pfcr(*12+*23)) 
" ''2p,3p 



(8) 



For the two pulse echo one should set i23 = in the above expression. 
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FIG. 3: Experiment (symbols) vs. theory (lines) for the time dependence of the two pulse echo amplitude in BK7. Straight 
lines (colored online) from the top to the bottom corresponds to experimental temperatures in ascending order, 7.5mK, 25mK, 
55mK and lOOmK, respectively. 



Consider the qualitative nature of the echo amplitude time dependence, Eq. 



at various observation times. At 



short times, tj^ or ii2ii3 < ksTPoUokT 
order of magnitude during times ti2 ~ 



all TLS contribute with q-dependent weight to the signal, decaying by an 

h 



ksTPoUokT 



(two pulse echo) or ^23 



ksTPaUokTti 



(three pulse echo) in a 



qualitative agreement with earlier observations (see Reviews^ and references therein). At longer times the echo signal 



is determined by almost fully symmetric two level systems A ^ Aq ~ i?o, meaning g ^ 1 in Eq. 



Then one can 
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set 1 — 1 and expand the integral in Eq. ([8]) over q to infinity. This yields 

-^2p,3p oc ^- — — "2 . (9) 

ti2 In (1 + r/2p,3p(tl2 + hz)) (I + kres{tl2 + ^23)) 

Here one can determine two distinguishable behaviors (cf. Fig. |4|) depending whether the characteristic times 
are shorter (intermediate regime) or longer (long time limit) than the thermal or resonant TLS relaxation times 
l/kx, 1/kres- In the intermediate regime the echo signal is due to typical symmetric resonant TLS, q <^ I, \'y\ ~ jq. 
The power law time dependent echo decay is defined by the phase volume of these TLS < i which 

results in the laws (^12^13)"^ for two and three pulse echoes, respectively. At longer times the echo decay is 

determined by the relaxation and decoherence free subset of TLS weakly coupled to phonons, Tr(7^)|7o|~^ < jt^t^j- 

Their phase volume decreases with the time as t^g^ which determines the law of echo decay close to the experimentally 
observed one. Fig. |4l One should notice that three pulse echo data cannot be fitted by the exponential decay 
occurring if the distribution of the TLS-phonons coupling constant, Eq. is ignored. The crossover between the 
two quasi-power law regimes is used below to estimate the TLS relaxation rates. 

One should notice that the "longitudinal" p art of TLS elastic te nsor cannot be ignored in the very long time limit, 
^23 ^ 1/fcros- It leads to an extra factor l/-\/5/2 + fc;(ti2 +^23), where ki ~ k^cs{vt/vi)^ ^ O.lfcics for the phonon 
stimulated relaxation and changes the long time asymptotic behavior from / oc t^^ to t^^'^ at ti3 > 10/fcios- Since 
in the experiment under consideration the times of measurement are shorter (see Fig. |4] and the relaxation time 
estimates) the longitudinal part of elastic tensor is not significant. 
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FIG. 4: Experiment (symbols) vs. theory (lines) for the time dependence of the three pulse echo signal in BK7. Straight 
lines (colored on line) from the top to the bottom corresponds to experimental temperatures and times, t\2, in ascending order, 
15mK (2/xs, 5.5/is, 15/is), 40mK (2/is, 5.5/is, 15/^s), respectively. 

Theoretical model is expressed by Eq. (jS]) and experimental data are available for two and three pulse echoes in 
Borosilicate glass BK7-i2ii^ for the wide range of temperatures including T = 7.5, 15, 25, 40, 55, 70 and lOOmK for 
two pulse echoes (l/is< ti2 < 1ms) and T = 7.5, 15, 25, 40, 55, 70mK and durations ti2 = 2, 5.5 and 15/iS, for three 
pulse echoes (6/xs< t23 ^ 10ms). Applying simple Monte-Carlo algorithm minimizing the mean squared deviations 
between logarithms of experimental data and theoretical fits to find an optimum set of the free parameters defined 
above (see Appendix for detail), we obtained a good fit of all experimental data choosing A — 4.66 • IG^s^^K^'^, 
B = 6.3 • 10*s~^K~^, a = 0.018, ^ = 0.46, 772P = 0.63, ?73p = 9.45. Some examples of comparison of theory vs. 
experiment for two-pulse and three pulse echo time dependencies are shown in Figs. [3] and S) For all data sets the 
accuracy of the fit is close to the accuracy of the experiment (see Appendix) . 

Let us compare our fitting parameters with available experimental data for BK7. The latter three parameters 
^ = 0.46, 7y2p = 0.67, rysp = 9.45, are of order of unity as expected. Difference between parameters ri2p = 0.63 and 
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Vsp = 9.45 can be because of the different scenarios of decoherence for two and three pulse echoes (see Appendix). 
The more accurate study of the problem requires the application of the telegraph process formalism-- which is outside 
of the scope of the present paper. 

The first three parameters, A, B and a can be directly compared with the results of the previous investigations 
of BK7. For convenience we assign the subscript "echo" to our Monte-Carlo estimates, while the estimates based on 
previous studies are given with the subscript "p" . According to the classical paper— one can express the parameter 
Ap from Eq. ([1]), in terms of the TLS coupling strengths to longitudinal, 7;, and transverse, 74, phonons having 
velocities ci and Cj, respectively, as 



Using the earlier experimental results for BK7— we get Ap 4.0 • lO^s^^K^'^, which is consistent with the estimate 
based on the two and three pulse echo experimental data, Aecho = 4.66 • lO^s^^K^^. Some overestimate of the 
relaxation rate can be due to the small effect of longitudinal part of TLS-phonon interaction and/or dipole-dipole 
TLS interaction, as noticed above. 

Consider the relative strength of TLS-TLS interaction, x = PoUq. It is defined by the fitting parameter accho = 
— 0.018 as Xccho ~ 0.44 • 10~^. Alternatively, one can estimate the average TLS-TLS interaction constant Uop 
using the TLS coupling with phonons, jp determined by the relaxation parameter Ap. Applying the Fermi golden 

rule to the symmetric TLS relaxation we get Ap — (^j^ + ^'s'^jnAp ■ Comparing this expression with the above 

estimates of the parameter Aodio.p one can find the average elastic tensor squared invariant as?^ |7p| = 4.18eV. TLS 
coupling constant, Uq^ is defined as a quadratic form of two TLS elastics tensors depending on their orientations. 
Assuming 7's to be traceless random Gaussian matrices with I7I obeying the statistics, Eq. (0), we evaluated Uq 
numerically and obtained Xp ~ 0.5 • 10^'^ consistently with the echo based estimate. 

Consider the TLS-TLS interaction contribution to the TLS relaxation rate characterized by the parameter B 
(see Eq. (U])). Using the previously estimated diniensionless parameter x a-nd the qualitative expression of Rcf.^-, 
Bth = lOfcsX^/^K^^s^^, we get Bth ~ 164K^^s^^, which is smaller than the echo based estimate, Bccho ^ 64000, 
by the factor of 380. To attain an agreement between experiment and theory one should set Xth ~ 0.0035 which is a 
factor of 7 larger than the estimate above, x ~ 0.0005. 

What is the source of this discrepancy? One possible explanation is that there is a different TLS relaxation 
mechanism which leads to the faster relaxation than the one proposed in Refs.^**'^^ Another reason could be that the 
qualitative theoryi^ misses a large numerical parameter, which can account for the above discrepancy. This looks 
quite realistic if we consider how the parameter x enters other theories. For instance in the expression for the spectral 
diffusion induced phase decoherence, Eq. ([7]), it enters together with the large numerical factor 7r^/24 ^ 40. This 
factor comes from the unity sphere volume contribution and integration over energies. One can expect a similar rise of 
the effective parameter x in the theory of interaction stimulated relaxation, which is more than sufficient to interpret 
the present experimental data. Of course, it is desirable to develop a more rigorous theory of the TLS-TLS interaction 
stimulated relaxation, which is beyond the scope of the present paper. 

Alternatively the theory can be tested by the experimental verification of the relaxation rate parametric depen- 
dence k oc x^- For instance, using Ref;^ one can estimate x ^ 0.00042 for a— Si02. Accordingly the interaction 
stimulated relaxation rate in a— Si02 should be smaller than that in BK7 by a factor of 2, which is interesting to 
verify experimentally. Preliminary analysis of the short time two-pulse echo measurements performed using the same 
experimental setupi^ yields -Bsi02 ~ 22000 w i?BK7/3 which is approximately consistent with our expectation. 

Thus we proposed a general theory interpreting two and three pulse dielectric echo measurements at microwave 
frequency, wq, in low temperature amorphous solids. We show that at long times the echo amplitude is determined 
by the relaxation and decoherence free subset of symmetric TLS having small coupling with phonons, Tr(7^) ^ 1/^23, 
with the phase volume approximately scaling as 1/^23- This power law describes the decay of the echo with the time in 
the long time limit. Our theory interprets the recent experimental data in BK7 for two and three pulse echoes within 
the experimental accuracy. One can conclude that the dielectric echo measurements can be used to determine the 
TLS relaxation rate based on the crossover between different power law time dependencies rather than the exponential 
decays. 

The relaxation rate temperature dependence satisfying the experimental data can be represented as a superposition 
of two contributions due to phonons, kph oc T^, and due to TLS-TLS interaction, kxLS T. This observation is 
consistent with previous measurements. The relaxation rate due to TLS - phonon interaction extracted from the echo 
experimental data is consistent with the previous measurements in BK7. The nature of the quantitative discrepancy 
between experiment and theory for TLS-TLS interaction stimulated relaxation in glasses is discussed and further 
experimental verifications of that theory are proposed. 
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I. APPENDIX 



A. Echo decay due to the spectral diffusion 



Here we characterize the echo decay associated with the TLS phase fluctuations caused by its interaction with 
neighboring TLS3^ This decay is defined by the following average function (for two-pulse echo one should set t23 — ^12) 



exp 



dT{E{T)-E{T+t23)) 



(11) 



where E{t) is the time-dependent energy of the TLS under consideration having average asymmetry A and elastic 
interaction tensor 7. If a TLS energy is time independent, then two phases in the exponent cancel each other and no 
decay occurs. The time dependent part of TLS energy can be represented asi^i^ 



6Eit)=qxJ2u,S^it), 



E ' 



70 



(12) 



where qxUj is the interaction of the given TLS with some TLS j, factors q and x are separated to have the bare 
interaction Uj independent of the TLS under consideration similarly to Refs.^^. Pseudospin 1/2 operators S^{t) 
describes the present states of neighboring TLS being 1/2 in the ground state and —1/2 in the excited state. The 
spontaneous transitions between these states determine the time dependence of TLS energy. 

Then ignoring correlations between weakly interacting TLS one can represent the phase factor, Eq. (|lip . as 



exp 



^ Jo 



dT{E{T)-E{T + h3)) 



J|exp 



/ ixqUj 



(13) 



Using the general method of averaging the products over neighbors^^'^^ characterized by the density of states 



R 



exp 



<l[F{i)>=Yl{l-dR dEi dAoiP(£;i,Aoi)(l-F(£;i,Aoi,R)) 



dR dEi dAoiP(Si,Aoi)(l-F(^i,Aoi,R)) 



one can reexpress Eq. ([T^ as exp [~qxG{t 12, t23)] where the function G is defined as 



G(ti2,i23) = X dx 



dA, 



01 



Aoi 



+ 00 



An 



dEi I 1 — ( exp 



2hx3 



dT{S!{T)-St{T + h2)) 

X = PoUa. 



(14) 



(15) 



Here the replacement of variable x^ = has been used. This substitution results in the prefactor of the average 

absolute value of TLS coupling constant, Uq. 

In the case ti2 <C ^23, realized in the majority of three pulse echo measurements one can simplify the integral in 
the exponent in Eq. as jj" dT{SS^{T) — SS^{T + t23)) « ii2(5'^(0) — 5'^(i23))- The exponent can be then evaluated 
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explicitly using the properties of the spin 1/2 operator, = cos(a/2) + 2i5'^ sin(a/2), and the spin-spin correlation 
function < 6S^{t)SS^{0) TtV^' "^^ere k = k,{E,T) is a TLS relaxation rate and k,{E,T) is a 

4 cosh 2kgT ) 

relaxation rate at temperature T for symmetric TLS (A = 0) with energy E. Using these definitions one can express 
Eq. ([T5|) in the form 

^ dv /■+°° dEi -v^k4Ei,T)t23 

w = — . (16) 



The integral over x can be evaluated exactly as / rfx (l — cos (^)) — ^f-|ii2|- Then, after the replacement of variable. 



3?i 7o Jo cosh (li) 



The analytical evaluation of the integrals in Eq. (|T6|) is not possible. There are two asymptotic behaviors for long 
and short delay times i23 compared to the thermal TLS relaxation time Tit ^ ks{kBT,T)~^ , namely. 

. ^kTlt23, ksikBT,T)t23 < I, 

^{tl2,l23) ~ I U2UT^^^ks{kBT,T)t23), k,{kBT,T)t23 > 1, 

Ur^^-^kBTP^U., /^- ^.(2.feg,T).. _ ^^^^ 

o Jo cosh [u) 

According to the earlier studiesi^"— TLS relaxation rate in dielectric glasses can be approximately represented as the 
sum of two contributions 

eV , f E 



associated with the TLS-phonon and TLS-TLS interactions, respectively. The phonon contribution to fc^i can be 
evaluated analytically22i as kxiph = ^AT^, while we represent the average TLS-TLS interaction contribution in the 

form kriint — '^S.AT'^ with some unknown fitting parameter of the theory, ^ '-^ 1, depending on the specific of the 
interaction stimulated relaxation. 

For the analysis of experimental data we used simplified analytical expressions 

akBTti2\n{l + ri2pkTti2) 

Or2piti2j = 7 ; 

n ?72p 

akBTti2 In (1 -t- rispkritu + %)) _ 

kr - AT^ + iBT. a = ^Pof/c (20) 

which handles all asymptotic behaviors for the three pulse echo, Eq. (ITS)) , and should be valid at least qualitatively 
for the two pulse echo (cf. Ref.— ). The time ^23 in Eq. ([TS)) is replaced with the sum ti2 + <23, which can be justified 
by the exact perturbation theorj^ at intermediate times t23 < k^^. To emphasize the difference between two cases 
we introduced two different free parameters ?72p and 773p in two and three pulse echo expressions for the decay rates in 
Eq. ([T])- Unknown parameters a, ??2.3p, A^ B, ^ are all used as free parameters of the theory determined by the best 
fit of the experimental data. Here we show all experimental data skipped in the main body of the paper together with 
theoretical curves for two pulse (Fig. [5]) and three pulse (Fig. ^ echo studies at different temperatures to illustrate 
nearly perfect agreement between the experiment and the theory. 



G3p{tl2,t23) 



B. Monte-Carlo search for the free parameters 



As an input we have the sequence of experimental data yi for two and three pulse echoes in a wide range of 
temperatures and times of measurements to be fitted by the theoretical expression, Zi(P) (Eq. (8) in the main text). 
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FIG. 5: Experiment (symbols) vs. theory (lines) for the time dependence of the two pulse echo amplitudes in BK7. Straight 
Unes (colored online) from the top to the bottom corresponds to experimental temperatures in ascending order, 7.5mK, 15mK, 
40mK and 70mK, respectively. 



depending on the set of unknown fitting parameters P. 
root mean squared deviation of experiment and theory 



In each step of the procedure we calculate the logarithmic 



(21) 



The scaling constant c is determined by the minimum requirement for the function Eq. (j2ip . The logarithmic deviation 
criterion has been chosen because experimental echo amplitudes change several orders of magnitude depending on the 
time of measurements. 

Then the new set of parameters, P', is created multiplying each parameter by the random factor R — e^^ with a 
random number x distributed uniformly between —1 and 1. Initially we choose ijj = 2. For this new set of parameters 
the new deviation D' is calculated using Eq. ((2T|) . In case oi D' < D we took P' = P, or left P without changes 
otherwise. The procedure is repeated 1000 times. Then we reduce ip by the factor of 2 and repeated the search again. 
This procedure is continued until ip — 0.125 and the result is considered as a final answer for the set of optimum 
parameters. We applied this algorithm several times obtaining the final answers different by no more than 10% from 
each other. 

The application of this method yields an accurate fit of all experimental data. Most strong deviations occurs for 
the two-pulse echo at short times ti2 ~ 1/iS. Perhaps this is because the driving field amplitude is not so small 
and perturbation theory used to derive Eq. (8) in the main text is not quite applicable. This is not important for 
intermediate and long times where only "symmetric" TLS contribute to echo, where the pulses related part of echo 
amplitude is decoupled from its decay. 
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(colored online) from the top to the bottom corresponds to experimental temperatures in ascending order and different times 
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